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Abstract 



The phase diagram of the Blume-Capel model on a semi-infinite simple cubic 



lattice with a (100) free surface is studied in the pair approximation of the 
cluster variation method. Six main topologies are found, of which two are 



<^> new, due to the occurrence of a first order surface transition in the phase 

| with ordered bulk, separating two phases with large and small surface order 

| parameters. The latter is a new phase and is studied in some detail, giving the 

behaviour of the order parameter profiles in two typical cases. A comparison 
is made with the results of a low temperature expansion, where these are 
q | available, showing a great increase in accuracy with respect to the mean field 



approximation . 
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I. INTRODUCTION 



The Blume-Capel (BC) model []l],|2| is a spin-1 Ising model with single-ion anisotropy 
that has been introduced as a model for magnetic systems and then applied to multicompo- 
nent fluids 0] and extended to the Blume-Emery-Grifliths model 0] for He 3 -He 4 mixtures. 
It is not exactly solvable in more than one dimension, but it is has been studied over infinite 
rf-dimensional lattices by means of many different approximate techniques, and its phase 
diagram is well known. 

In recent years, when the theory of surface critical phenomena started developing, some 
attention has been devoted to the study of the BC model over semi-infinite lattices, with 
modified surface couplings. In particular Benyoussef, Boccara and Saber || have determined 
the phase diagram in the mean field approximation, reporting four possible topologies at 
fixed bulk/surface coupling ratios, while Benyoussef, Boccara and el Bouziani |J have done 
a similar analysis using a real space renormalization group transformation. Other works, 
referring to particular regions of the phase space, are those by Peliti and Leibler |7[ and 
Crisanti and Peliti || (real space renormalization group), Jiang and Giri || (mean field 
approximation), Tamura JT(| (effective field approximation) and the present authors |TT|JT2] 
(CVM, low temperature expansion). 

All these results show that, as in the case of the spin-1/2 Ising model JIB] , it is possible 
to have a phase with ordered surface and disordered bulk, which is separated from the 
completely ordered phase by the so-called extraordinary transition and from the completely 
disordered phase by the surface transition. When such a phase is absent, the transition 
between the completely ordered and the completely disordered phase is named ordinary and 
the meeting point of the lines of these three kinds of phase transitions is said special, and 
is generally a multicritical point. While it is well established that the ordinary transition 
can be either first or second order, the possibility of a first order surface or extraordinary 
transition has been ruled out by Benyoussef, Boccara and el Bouziani || and Tamura ||10|| . 
but not by the other authors. 
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In the present paper, which can be regarded as a refinement of a previous work |TT 
we apply the pair approximation of the cluster variation method (CVM) fll4]-|T6|j to the 
Blume-Capel model defined on a semi-infinite simple cubic lattice with a (100) free surface. 
The dependence of the local quantities (order parameters, energy density, and so on) on the 
distance from the free surface is taken into account by dividing the lattice into an infinite 
set of layers parallel to the surface and then treating the Nth layer (we will use N = 5) as 
if it was bulk. 

Our results confirm the existence of first order extraordinary and surface transitions (for 
the surface one we obtain very good agreement with a previous low temperature expansion 
analysis). Moreover, we show clearly that another first order transition can occur at the 
surface in presence of an ordered bulk, which separates two phases with large and small values 
of the surface order parameters, the latter being a new phase. When this new transition is 
included in the phase diagram, the classification given by Benyoussef, Boccara and Saber 
H is completed with the addition of two new topologies. 

The paper is organized as follows: in section II we define the model and determine its 
ground state; in section III we give the pair CVM free energy functional and describe the 
nested iteration procedure we use to find its minimum; in section IV we describe our results, 
giving the complete classification of the phase diagram topologies and the behaviour of 
the order parameters and their profiles in the new phase. Finally, in section V, we briefly 
summarize and discuss our results. 

II. THE MODEL AND ITS GROUND STATE 



The BC model on a semi-infinite lattice A with a free surface denoted by <9A is defined 
by the (reduced) hamiltonian 

(3H = -J £ 8 iSj -J B £ S ^ + A$> 4 2 + A s $> 2 fc , (1) 

(ij}<£dA (kl)cdA igdA kedA 

where i,j,k,l are site labels, Sj = ±1,0 is the ^-component of a spin-1 operator at site 
i, J and J s denote the bulk and surface nearest neighbors (reduced) interactions (assumed 



positive) , respectively, A and A s denote the bulk and surface single ion (reduced) anisotropy 
and (3 = (k-gT)" 1 , with k-Q Boltzmann's constant and T absolute temperature (from now 
on, we will express the temperature in units of the bulk interaction, by setting T = 1/ J) . 
Furthermore, we introduce the coupling ratios R— J/J s and D = A/A s . 

As it is often the case with classical discrete models, we can easily determine the ground 
state of our model by looking for the lowest energy configurations. This analysis is very 
important since it drives all subsequent considerations. 

For the bulk, as is well known, we have a state with broken spin-flip symmetry for 
A < 3 J, given by e.g. Sj = +1 for all bulk spins, and a state with unbroken symmetry, 
Si = 0, for A > 3 J. These two states are associated with ordered and disordered bulk 
phases, respectively, at finite temperature. 

At the surface the situation is slightly more complex. Let us define the states LO and 
HO by Si = 0,Wi G <9A and Sj = +l,Vz G OA, respectively. The corresponding energies per 
site will be E LO = and E^o = A s — 2J S — m J, where m = in the bulk disordered 
phase and 1 in the bulk ordered phase (m is but the bulk order parameter (sj) at zero 
temperature). 

We can now distinguish the following cases JT7J: 



a) when D < 3/2/(2 + R) there is a surface transition at A/ J = <5 = D(2 + R)/R, 
in presence of an ordered bulk, and the states LO and HO are associated, at finite 
temperature, to two phases with ordered surface, but with different values of the 
surface order parameters: the phase HO (we use the same symbols for the phases and 
the corresponding pure states at T = 0) has large order parameters and is the usual 
ordered phase, while LO is a new phase with small order parameters; this situation has 
not been considered before, except in the limit A — ► — oo, which yields the spin-1/2 
Ising model with a spin-1 free surface studied by Kaneyoshi |T% , in which a first order 



surface transition in the presence of ordered bulk is obtained, but neither studied in 
detail nor reported in the phase diagram; 
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b) when 3R/(2 + R) < D < 3R/2 there is a surface transition at A = 3J, in correspon- 

dence of the bulk transition, and the states LO and HO are associated to phases with 
disordered and ordered surface, respectively; 

c) when D > 3R/2 there is a surface transition at A/ J = 5\ = 2D/R, in presence of a 

disordered bulk, and the states LO and HO have the same meaning as in the previous 
case. 



III. THE CVM PAIR APPROXIMATION FOR SEMI INFINITE SYSTEMS 

In order to study the finite temperature properties of the model, we introduce now the 
pair approximation of the cluster variation method. 

The CVM is based on an approximate expression of the entropy of the model as a sum 
of contributions by all the elements of a set M. consisting of certain maximal clusters and 
all their subclusters. In the pair approximation, one takes as maximal clusters all nearest 
neighbor pairs, and thus the only subclusters to be considered are the sites. To each cluster 
7 is associated a cluster (reduced) entropy given by 

S 7 = -Tr(p 7 lnp 7 ), (2) 

where p 7 is a cluster density matrix, to be determined minimizing the free energy. The total 
reduced entropy is then approximated by 

S — ^ Cl-yiy^ S '-y , 



where the coefficients a 7 obey [|i~5l . |TG 



a 7 = 1, \/aeM (4) 

and z/ 7 represents the multiplicity (number of clusters per site) of the cluster 7. 

Since we have to deal with local quantities (density matrix elements, in our formulation) 
which depend on the distance from the surface, we divide our lattice, as mentioned in the 



introduction, into two dimensional layers parallel to the surface, labeled by an integer n, 
n — 1 being the surface layer. Furthermore, to avoid working with an infinite number of 
variables, we will treat the iVth layer (typically we will choose iV = 5) as if it was bulk, 
by means of suitable constraints, and the bulk itself will be considered in the CVM pair 
approximation, too, using the results obtained in [ITS . 



Following |]17| ) we introduce the density matrices p 2 n \si, s 2 ) for the pairs having both 
sites in the nth layer, p 2 ™' n+1 ^(si, s 2) for those pairs having spin s\ in layer n and spin s 2 
in layer n + 1, and Pi (si) for sites in layer n. The pair density matrices must satisfy the 
obvious constraints 



S2 S2 S2 

which ensure that the site density matrices can be properly defined, e.g. by 

^\s 1 ) = W[p^(s 1 ,s 2 )+p^(s 2 ,s 1 )}. (6) 

Finally, in order to make the iVth layer represent the bulk in an effective way, we will treat 
the corresponding pair density matrix as a constant, defined by 

P { 2 N \si,s 2 ) = Pbulk{si,S 2 ), (7) 

where Pbuik is the normalized bulk pair density matrix given by the CVM pair approximation 
|T9| . Notice that the normalization of pbuik and the constraints Eq. |5| imply the normalization 
of all the density matrices we use. 



Using the coefficients a 7 and z/ 7 reported in [|IT[] one readily obtains the (reduced) free 
energy density 



/ = -2J s Tr( SlS2 p«) + A s Tr( S ^ 1} ) - J T Fr(s 1 s 2 ^ ) 

N-l 

+ [-2JTr( Sl s 2 p ( 2 n) ) + ATr(^pW) _ JTr( Sl s 2 p^ n+lh 



71=2 

N-l 



+ £ [2Tr(pf lnp 2 n) ) + Tr(pf'" +1) \np^ n+l) ) + a^Ti(p^ lnp^ 

n=l 
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N-l 

si n=l 
iV 

si n=2 



P2 



(n '" +1) (Sl,S 2 )-pi n) (^ 



Ep^ 1,B Wi)-rf°(«i) 



*2 



(8) 



where A± are Lagrange multipliers, ajj 1 = —4 and af^ = —5 for n > 1. 

Minimizing this functional corresponds to solving the following equations, in the form of 
the natural iteration method 12T 



P2 {si,s 2 ) = exp 



x exp 



(A^O + A^)) 

A, 
4 



x 



; — s / 2 , 2> 
JsSlS 2 + S 2J 



(n>l)/ x 

p 2 '{s 1 ,s 2 ) = exp 



X 



x exp 



A 



Jsi-s 2 - —(s 1 + s 2 ) 



5/4 



(n,n+l)/ \ 

p 2 ; (si,s 2 )=exp 



A H (S1 )_ A ^ (S2 ) exp(J SlS2 ), 



fn+l) 



(9) 



together with the minor iteration equations [21,221 



E(n,n+1)/ \ 
p 2 ; (si,s 2 ) 



Ai n) ( S i) 



r+l 



1 



Pi (si 



A W ( Sl ) 



r+l 



AL n) ( Sl ) 



Ell 
p 2 



(n— l,n) 



S2,Slj 



+ -ln^- 

r 2 



(n) / \ 
Pi (Si/ 



(10) 



for the Lagrange multipliers, where r and r + l denote successive steps of the iteration 
procedure. 

The method of solution can be described as follows: 

0. initialize all pair density matrices with a rough estimate of the solution; 

1. determine the site density matrices according to Eq. ||; 

2. recalculate the pair density matrices according to Eq. ||with all Lagrange multipliers set 

to 0; 



3. solve Eq. [10| for the Lagrange multipliers by simple iteration; 



4. determine the pair density matrices using Eq. ^| with the values of the Lagrange multi- 

pliers obtained in step 3; 

5. iterate steps 1-4 until the desired precision is reached. 

Recalling that, in region of the phase space close enough to first order transition lines, 
our free energy should have different minima corresponding to stable and metastable phases, 
one has to repeat the above procedure with different initializations, corresponding to the 
different phase, and then choose the solution of minimum free energy. Once the solution is 
found one readily obtains the order parameters 

™»= E sPi\s) (11) 

s=±l,0 

Qn= E * 2 Pi n) (*) (12) 
s=±l,0 

and, in an analogous way, all the correlation functions. 



IV. PHASE DIAGRAM 

By means of the nested iteration procedure described in the previous section, we have 
determined the phase diagram of our model in the plane (A, T = 1/J), for fixed values of 
R and D. We have obtained six main topologies, reported in Fig. 1, extending thereby 
the classification by Benyoussef, Boccara and Saber H], in which cases (a) and (d) were 
missing. These topologies are identified by specifying the ground state, according to the 
classification given in section II, and the result in the Ising limit A, A s — > — oo. Recalling 
the known results about the critical behaviour of the semi-infinite Ising model |T3| 



one 



realizes that in this limit there are two possibilities: for R > R c the model exhibits an 
ordinary transition at the Ising bulk critical temperature, while for R < R c it undergoes an 
extraordinary transition at the bulk critical temperature and a surface transition at a higher 
temperature. The best estimate for R c is about 0.66 J23|,|4]]. In Fig. 1, case (a) corresponds 
to R > R c and D < 3R/(2 + R), case (b) to R > R c and 3R/(2 + R) < D < 3R/2, case (c) 
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to R > R c and D > 3R/2, case (d) to R < R c and D < 3R/{2 + R), case (e) to R < R c and 
3R/(2 + R) <D < 3R/2 and, finally, case (f) to R < R c and D > 3R/2. 

Notice that almost all of the above topologies, with the exception of cases (b) and (f), 
could be further divided into subtopologies based on the order of surface and bulk transitions 
at the special point M and on the eventual coalescence of that point with the critical point 
C. This is exemplified in Fig. 2, which refers to cases (a) and (d) of Fig. 1. 

It is also interesting to look at the behaviour of the order parameters in the cases (a) 
and (d) of Fig. 1, where the new LO phase appears. In Figs. 3a and 3b, which refer to case 
(a), we see that the surface layer order parameters increase from zero at low temperatures 
(it can be shown exactly that at low temperatures they behave as exp(J — A s )), while the 
order parameters of the inner layers behave nearly as the bulk ones. Furthermore, in Figs. 
4a and 4b, the effect of the new surface phase transition is clearly shown. 

Let us now turn our attention to the first order surface transition which occurs in pres- 
ence of a disordered bulk, which existence, questioned by previous works [PJlQfl, is clearly 
confirmed in our approximation (Fig. 1, cases (c) and (f)). In Fig. 5 we have compared, 
for one choice of R and D, our results with those obtained by mean field approximation || 
and low temperature expansion ||12|| . It is clearly seen that our approximation is a marked 
improvement with respect to the simple mean field approximation, and the very good agree- 
ment with the low temperature expansion results is a strong argument in favour of the 
existence of such a transition. 

Finally, we have not found any reentrant phenomenon in the phase diagram, contrarily 
to previous results This discrepancy should be due to the very rough approximation 



used in |[L 1|| , where only one layer above the bulk was considered and, a factor even more 
important, the method of solution adopted did not satisfy the first of the constraints Eq. |5| 
in the case n = N (=2, in that case). 
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V. CONCLUSIONS 



We have determined the phase diagram of the semi-infinite Blume-Capel model using 
the pair approximation of the cluster variation method. Our analysis has shown, at fixed 
R and D, six main topologies, two of which are completely new because of the appearance 
of a new ordered phase and, consequently, of a new surface first order phase transition in 
presence of an ordered bulk. This new phase has been studied in some detail and typical 
behaviours of the order parameters profiles have been reported. We have also given an 
example of how these main topologies may be slightly modified by varying R and D. 

A question concerning the existence of the surface first order transition in the presence of 
disordered bulk has been addressed and clarified by comparing our results with previous re- 
sults obtained by a low temperature expansion and, finally, the issue of reentrant phenomena 
in the surface transitions, raised in a previous paper, has been shortly discussed. 
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FIGURE CAPTIONS 



FIG. 1. Phase diagram topologies: heavy and light lines denote bulk and surface transitions, 
respectively, while solid and dashed lines denote second and first order transitions; C and 
M denote critical and multicritical points. 

FIG. 2. Phase diagram for D = 0.2 and R = 1.00, 0.50, 0.32 and 0.25. 

FIG. 3a. Order parameters m n for R = 1, D = 1/2, A/J = 1.6; n — 1 (lowest curve), 2, 3, 
4 and 5 (bulk, highest curve). 

FIG. 3b. Order parameters q n in the case of Fig. 3a; n = 1 (lowest curve), 2, 3, 4 and 5 
(bulk, highest curve). 

FIG. 4a. The same as Fig. 3a for R = 1/2, D = 1/2, A/ J = 2.6. 
FIG. 4b. The same as Fig. 3b for R = 1/2, D = 1/2, A/J = 2.6. 

FIG. 5. First order surface transition line for it! = 1, D — 1.53 as given by present method 
(solid line), mean field approximation (long dashes) and low temperature expansion (short 
dashes) . 
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